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Abstract 


Jr 


t 


Conditions  are  given  for  uniqueness  of  limit  cycles  for  autonomous 
equations  in  the  plane.  The  results  are  applicable  to  codimension  two 
bifurcations  near  equilibrium  points  for  vector  fields. 
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SI.  INTRODUCTION 

Consider  a  system  of  autonomous  differential  equations  in  IR^  , 
d  ^  2,  or  in  some  infinite  dimensional  space.  Suppose  that  there  exists 
an  equilibrium  point  which  is  non-hyperbolic  and  is  doubly  degenerate. 

In  other  words,  it  is  a  codimensional  two  singularity.  More  specifically, 
we  assume  that  the  linear  variational  equation  in  the  center  manifold 
has  one  of  the  following  linear  parts 
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For  more  detail,  see  Arnold  [1,2],  Carr  [8],  Chow  and  Hale  [9] 
Guckenheimer  and  Holmes  [11]. 

•  The  classification  of  vector  fields  near  a  singularity  of 
r 

|  for  generic  perturbations  can  be  found  in  Arnold  [1],  Bogdanov 

: 

*  and,  for  symmetric  perturbations,  in  Carr  [8]  and  Takens  [14]. 
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classification  is  not  simple.  The  main  difficulty  is  to  prove 
uniqueness  of  the  limit  cycle.  It  is  usually  proved  by  using  proper¬ 
ties  of  elliptic  integrals  depending  on  a  parameter  (see  Arnold  [1], 
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Bogdanov  [1,2,3],  Carr  [8],  Chow  and  Hale  [9],  and  Guckenheimer  and 
Holmes  [11])  and  is  related  to  the  weakened  Hilbert  16th  problem 
(Arnold  [1],  p.  303).  In  [12,13],  Il'jasenko  showed  that  these 
properties  of  elliptic  integrals  could  also  be  obtained  by  methods  in 
algebraic  geometry. 

The  classification  of  vector  fields  near  a  singularity  of 
type  or  A^  is  far  from  being  complete.  Since  the  dimension  of  the 
system  is  >3,  many  new  types  of  dynamical  behavior  can  occur.  The 
first  step  in  any  attempt  at  classification  is  to  put  the  equations 
in  normal  form  and  then  try  to  classify  the  truncated  equations;  namely, 
the  normal  form  equations  up  through  polynomials  of  degree  <_  k  for 
some  fixed  integer  k.  Because  of  the  nature  of  the  linear  parts 
A2  or  A^,  these  equations  are  very  symmetrical  and  the  flow  can  be 
reduced  to  a  polynomial  equation  in  the  plane.  It  is  this  polynomial 
equation  which  is  the  subject  of  this  paper.  For  this  equation,  the 
existence  and  number  of  periodic  orbits  plays  a  fundamental  role  in 
the  classification. 

Most  of  the  paper  is  concerned  with  the  case  A^.  In  Section  2, 
for  the  case  A^  we  show  how  one  obtains  the  relevant  equation  in  the 
plane  and  point  out  that  the  bifurcation  diagram  is  easy  to  obtain 
everywhere  except  in  the  neighborhood  of  one  point  for  which  the 
unperturbed  equation  is  Hamiltonian.  The  discussion  of  period  orbits 
reduces  to  the  consideration  of  the  monotonicity  properties  of  a 
certain  function  which  is  the  ratio  of  Abelian  type  integrals.  For 
certain  values  of  the  parameter,  these  are  Abelian  integrals  and  one 
can  analyze  this  function  by  using  methods  from  the  theory  of  analytic 
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functions  and  algebraic  geometry.  In  Section  3,  we  discuss  a  case  of 
this  type  and  show  periodic  orbits  are  unique.  The  method  of  proof 
follows  earlier  ideas  used  by  Carr  [8].  Cushman  and  Sanders  [10]  also 
have  discussed  this  case  and  have  shown  very  clearly  the  connections 
with  the  Picard-Fuchs  theory  of  Abelian  integrals.  Bogdanov  presented 
this  result  in  [4]  with  a  different  proof.  Arnold  [1]  states  that  the 
method  of  Il'yashenko  [12,13]  can  be  applied  also  to  this  case. 

In  Section  4,  by  using  a  more  analytic  proof,  we  show  uniqueness 
in  case  for  situations  where  the  algebraic  methods  do  not  seem  to 
be  applicable. 

In  Section  5,  we  derive  the  bifurcation  function  for  case  A^ 
and  special  values  of  the  parameters.  The  general  case  is  still  open. 
Numerical  evidence  indicates  uniqueness  of  the  periodic  orbit. 

§2.  NORMAL  FORMS  AND  BIFURCATIONS 

In  this  section,  we  suppose  that  we  have  a  singularity  of  type  A^- 
In  order  to  classify  the  behavior  of  vector  fields  near  such  a  singu¬ 
larity,  consider  the  equation  in  R3 


u  =  Au  +  F (u) 


where  u  €  3R  , 


A  = 


-1 
,  0 


1  0 

a  0 

0  B 


2 

a,  B  €  R  are  small  parameters,  F  =  0(|u|  ),  as  j u |  0 ,  and  is 


smooth.  By  introducing  cylindrical  coordinates, 
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=  r  cos  0 
=  r  sin  0 

u3  =  5 

and  applying  the  theory  of  normal  forms  (method  of  averaging)  (see, 
for  example,  [9]),  one  may  assume  without  loss  of  generality  that 
our  differential  system  near  u  =  0  has  the  form 

9=1  +  fk(r2,s,ct,e)  +  0(k  +  l)/r 

r  =  rgk(r2's,“,3)  +  °(k  +  1 ) 

•  2 

s  =  h^fr  ,s,a,6)  +  0(k  +  1) 

where  fk(0,0,0)  =  gk(0,0,0)  =  hk(0,0,0)  =  3hk(0,0,0)/3s  =  0,  the 
functions  fk,rgk,hk  are  polynomials  in  r,  z,  a,  0,  of  degree 
_<  k  and  0(k  +  1)  denotes  terms  which  have  order  k  +  1  uniformly 
in  0. 

One  important  observation  about  these  equations  is  that  the  terms 
up  through  order  k  are  independent  of  0  and  that  the  equations 
truncated  up  through  order  k  always  have  the  manifold  r  =  0  in¬ 
variant  under  the  flow. 

An  attempt  to  analyze  the  above  equations  consists  of  the  follow¬ 
ing  steps.  First,  ignore  the  terms  of  order  k  +  1  and  analyze  the 
truncated  equations  completely.  Since  the  truncated  equations  are  in¬ 
dependent  of  0,  one  can  replace  the  time  variable  by  0  and  obtain 
an  equation  in  the  plane.  Because  the  truncated  equations  correspond 
to  a  planar  system,  periodic  orbits  play  a  major  role  in  the  discussion 
of  the  flow  for  a,0  parameters  varying  in  a  neighborhood  of  zero. 
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For  k  =  3,  we  will  give  below  a  complete  discussion  of  the  periodic 
orbits. 

The  next  step  in  the  analysis  is  to  discuss  the  effect  of  pertur¬ 
bation  terms  0(k  +  1).  These  terms  depend  upon  9  and  also  may  not 
possess  the  symmetry  properties  in  r  that  the  truncated  equations 
possess.  This  creates  significant  difficulties.  A  complete  solution 
to  the  complete  equations  is  not  available  at  this  time.  We  do  not 
discuss  this  problem,  but  mention  only  that  significant  progress  has 
been  made  by  Broer  [6],  Broer  and  Vetger  [7]. 

If  we  truncate  the  equations,  replace  t  by  0  and  keep  only 
terms  up  through  the  third  degree,  we  obtain  the  system 

—  —  2  —  2 
r'  =  r(a  +  as  +  dr  +  es  ) 

—  2  —  2  —  2  —  3 

s'  =  Bs  +  bs  +  cr  +  fr  s  +  gs 

where  "  '  "  =  d/d0  and  a,  b,  c,  d,  e,  f,  g  are  fixed  constants. 

The  tei'ms  neglected  from  the  truncated  equations  can  be  shown  to  be 
insignificant  (see  [9]). 

To  obtain  a  simpler  form  for  these  equations  we  rescale  the 
variables  by  letting 

r  =  Ibj^x 
s  =  IFlS^ 

9  =  -  zrirr  t 

b  |  e  |  ^ 

a  =  -b | c a 
e  =  -b  |  c  I"4  6 


to  obtain  the  following: 

.  -  3  ~  2 

x  =  ax  +  Bxy  +  dx  +  exy 

•  2  - 2-2  ~3 

y  =  By  -  y  -  (sgnbc)x  +  fx  y  +  gy 

where  B  =  -a/b  and  d,  e,  f,  g  are  constants. 

For  the  case  where  B  /  0,  b  c  <  0,  the  bifurcation  diagram  in  the 
(a, (3)  plane  is  easy  to  obtain  and  qualitatively  does  not  depend  upon 
the  cubic  terms.  The  bifurcations  consist  of  elementary  bifurcations 
of  equilibrium  points.  For  the  complete  equation,  this  means  only 
Hopf  bifurcations  and  the  coalescing  of  periodic  orbits.  For  a  complete 
discussion,  see  [9]. 

The  most  interesting  case  is  when  be  >  0.  It  is  convenient  to 

introduce  a  change  of  variables.  Choose  y  close  to  B/2  so  that 

B  -  2y  +  3d4Y  =  0  and  the  transformation  y  |+  y  +  y  together  with  a 

rescaling  of  x,y  leads  to  the  more  symmetric  form 

.  -  3  2 

x  =  ax  +  Bxy  +  d^x  +  d2xy 

*  B2  2  2,2  ,3 

y  =  —  -  y  -  x  +  d3x  y  +  d4y 


where  we  keep  the  same  labeling  for  the  constants,  even  though  they  are 
not  exactly  equal. 

If  we  perform  the  scalings 

x  ex,  y  -*■  ey,  6  e,  a  cX,  t  -+■  e  *t 
the  new  equations  become 

x  =  Xx  +  Bxy  +  e[d.x3  +  d2xy2] 

(2.1) 


2 

x  + 


e[d3x2y  +  d4y3] 


and  are  in  normal  form. 
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Equation  (2.1)  must  be  discussed  for  all  x  >_  0,  y  £  1R  ,  A  6  1R  , 
BCR,  and  e  in  a  neighborhood  of  e  =  0. 

In  the  following,  we  assume  B  >  0. 

Let  us  first  discuss  equation  (2.1)  for  e  =  0;  that  is,  the 


equation 


x  =  Ax  +  Bxy 


•  1  2  2 
y  =  4  -  y  -  x 


There  are  always  two  equilibrium  solutions  (0,±l/2).  The  equilibrium 

point  (0,1/2)  is  a  hyperbolic  stable  node  if  A  <  -B/2  and  a  saddle 

point  if  A  >  -B/2.  The  point  (0, -1/2)  is  a  saddle  point  if  A  <  B/2 

2  2 

and  a  hyperbolic  unstable  node  if  A  >  B/2.  If  A  <  B“/4,  there  are 

2  '’1/7 

other  equilibrium  points,  (±(1/4  -  A  /B")  ,  -  A/B) .  At  the  values 

A  =  ±B/2,  there  is  a  bifurcation  of  an  equilibrium  point  into  three 

equilibria,  each  of  which  is  hyperbolic.  This  bifurcation  still  remains 

2  2  l/'> 

when  e  /  0  is  small.  The  equilibrium  point  (±(1/4  -  A  /B  )  ,  -A/B) 

changes  its  stability  properties  at  A  =  0,  being  a  stable  focus  for 
A  >  0  and  an  unstable  focus  for  A  <  0.  The  point  A  =  0  thus  becomes 
another  place  where  bifurcations  can  occur. 

It  remains  to  analyze  the  behavior  of  the  solutions  of  (2.1)  for 
A  =  0,  e  =  0.  For  A  =  0,  e  =  0;  that  is,  the  equation 


(2.2) 


there  is  a  first  integral 


x  =  Bxy 


y  =  4  ~  y 


i 
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(2.3) 


H  =  £ 


[ 


A+l 


q+1  2 
-  xn  y 


xq+3 
1  +  B 


] 


where  1+q  =  2/B. 

The  fact  that  (2.1)  for  X  =  e  =  0  has  a  first  integral  implies 
that  cubic  terms  are  necessary  to  resolve  the  complete  bifurcation 
diagram  near  X  =  e  =  0.  To  obtain  this  bifurcation  diagram,  one  must 
discuss  the  periodic  orbits  of  (2.1)  for  (X,e)  small.  The  first 
integral  of  (2.2)  is  very  useful  in  such  an  analysis.  Let  us  briefly 
indicate  how  this  is  done. 

It  is  convenient  to  reparametrize  the  orbits  by  replacing  t  by 
xqt  to  obtain  the  system 

x  =  xq  [Xx  +  Bxy  +  efd^x3  +  d2xy2)] 

(2. 4) 

y  =  xq[i  -  x2  -  y2  +  e(d3x2y  +  d4y3)] 

The  derivative  of  H  along  the  solutions  of  (2.4)  satisfies 

dH  •  q+1.  • ,,  3  ,  q  • , ,  2  ,  3.  q 

(2 • 5)  fa  =  yxM  X+ey(d:x  +  d2xy  )xH  -  ex(d3x  y  +  d4y  )xM 


It  is  not  difficult  to  see  (Carr  [8  ])  that  a  necessary  and 
sufficient  condition  for  an  orbit  r  =  T(X,e)  of  (2.1)  to  be  periodic 
in  that 


(2.6) 


Using  (2.5)  and  (2.6),  one  obtains  a  bifurcation  function  for  periodic 
orbits  which  involves  Abelian  integrals.  (See  Sections  3,4  and 
Appendices  A,B).  It  is  then  shown  that  the  number  of  periodic  orbits 


is  then  reduced  to  the  discussion  of  the  monotonicity  properties  of  a 
scalar  valued  function  which  is  the  ratio  of  two  Abeli  *  integrals. 

Under  the  assumption  that  B  =  2,  d^  =  1 ,  d^  =  d^  =  d^  =  0  and  by 
assuming  the  monotonicity  of  the  above  function,  the  complete  bifur¬ 
cation  has  been  obtained  (see,  for  example,  Chow  and  Hale  [9]). 

In  the  following  section,  we  will  give  a  proof  of  this  monotonicity 
hypothesis.  Thus,  this  completes  the  bifurcation  diagram  for  equation 

(2.1)  with  B  =  2,  d^  =  1,  d2  =  d^  =  d4  =  0.  In  Section  4,  the  case 
B  >_  1/2  will  be  considered.  The  problem  of  monotonicity  is  still 
open  for  0  <  B  <  1/2.  In  Appendices  A,  B,  we  show  that  the  results 
obtained  are  still  valid  for  arbitrary  but  fixed  values  of  the  d^ 

not  lying  on  a  hyperplane  in  (dj ,d2 ,d3,d4)  space.  Note  that  this  is 
not  surprising  because,  on  a  Riemann  surface  of  finite  genus,  there  are 
only  a  fixed  number,  depending  on  the  genus,  of  linearly  independent 
holomorphic  differentials. 

§3.  UNIQUENESS  THEOREM  FOR  B  =  2 
Consider  the  equation 

x  =  Ax  +  2xy  +  ex 

(3.1) 


which  is  equation  (2.1)  with  B  =  2,  d1  =  1 ,  d2  =  d3  =  d4  =  0.  Let 


When  A  =  e  =  0,  (3.1)  is  a  Hamiltonian  system  with  the  Hamiltonian 
H.  Thus,  solutions  of  (3.1)  with  A  =  e  =  0  are  parametrized  by 
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where  c  £  1R  .  We  note  that  c  =  0  corresponds  to  the  heteroclinic 
orbit,  c  =  1/12  corresponds  to  the  fixed  point  (1/2,0)  while 
0  <  c  <  1/12  corresponds  to  a  periodic  orbit. 

As  remarked  in  the  previous  section,  a  necessary  and  sufficient 
condition  for  an  orbit  T  =  r(A,e)  of  (3.1)  to  be  periodic  is  that 

I  H  dt  =  0 


or 


0  = 


y(Ax  +  ex  )dt 


y(Ax  +  3ex  x)dt 


where  T  >  0  is  the  minimal  period  of  the  periodic  orbit  r(A,e).  By 
using  the  above  equation  and  letting  A  =  3ep,  one  obtains  a  bifurcation 
function  G(yi,e,c)  for  periodic  orbits  which  for  e  =  0  is  given  by 

G(p,0,c)  =  pJQ(c)  +  J2(c)  ,  0  <  c  <  Y2  ’ 


where 

(3.2) 


(3.3) 


C 

10M  =  | 


ydx 


J2(c)  = 


2  , 
x  ydx 


y  *  *(i -  t-  f) 


1/2 


and  0  <  a1  <  1/2  and  1/2  <  a2  <  are  the  two  real  positive  roots 

of  the  equation 


for  0  <  c  <  1/12.  Let 


!_*:_£=  0 
4  3  x 


(3.4) 


J2(c) 


0  <  c  <  1/12 


It  is  shown  in  [  9]  that,  if  P'(c)  /  0,  then  equation  (2.1)  has 
a  unique  asymptotically  stable  limit  cycle  for  appropriate  values  of 
(A,e)  i  (0,0).  Our  main  result  in  this  section  is  the  following. 


Theorem  5. 1  Let  P(c)  be  defined  as  above.  We  have  P'(c)  >  0 
for  0  <  c  <  1/12. 


prove  this,  let  c^  =  c ^  =  /a^  and 


(3.5) 


/  w  w  \ 

=(  4  '  C  "  3  j 


6  \  1/2 


0  <  w  < 


Define 


(3.6) 


I  =  f  w2nR 

n  Jc. 


(w)dw. 


By  setting  w  =  x  and  using  (3.2),  (3.3),  we  have 


2I2  -  3r  2Io  ■  Jo  -  0  *  c  ‘  u 


Hence , 
(3.7) 


Me) 

PW  =  W 


0<c<iT 


We  will  use  (3.7)  to  show  Theorem  3.1. 


Lemma  3.2 


2  .  2n 


I'  (C)  =  -  -  [  — 

}  2  I  R( 
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Proof.  From  (3.5),  2RdR/dc  =  -1. 


Lemma  3.3 


Proof. 


Lemma  3.4 


lim  P(c)  =  y 
1  4 

c-^ 


l'2M 

lim  P(c)  =  lim 
.  .  1  1  l0(c) 


’Of 


lim  P(c)  =  ~ 
c-*0  16 


Proof.  Integration. 


Lemma  3.5  Let  1(c)  be  the  column  vector  (I^jI^.I^).  Then 
T(c)  satisfies  the  differential  equation 

(3.8)  1(c)  =  A(c)  I '(c) 


where 


A(c) 


0 


1 

4 

c 


0 


1 

6 


i- 

64  4 


and 
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A  =  det  A(c)  =  f  (c2  -  )c 


Proof.  Let  us  first  prove  that  Ig(c)  =  3cIq(c)/2 


IJCc)/4 


the  relation 


2R  =  j  w  -  2w5 
dw  2 


it  follows  that,  as  differential  forms, 

R2 

Rdw  =  -jj-  dw 


dw  f  1  /  „  dR  1  \  w2  1 

■t  l' *t-  cJ 

] 


1  dw  1  2 

3  wdR  +  X  6  W 


Integration  by  parts  yields 


Vc>  ■ 


-  T'oW  - 11!'')  '  2‘VC) 


which  is  the  desired  relation.  Let  us  now  shown  that 


6IX  =  6cIJ  -  1'2 


By  the  definition  of  1^  and  by  integration  by  parts,  we  have 

h-f2  »2(t--t)x 

C1 

fC2  (  )L_  \  dw 

=  J  l  3  '  12  ;  R  ' 


From 


Hence , 


Finally,  we  show  that 

64  I2  =  6c  1'  -  3  IJ  +  48  cl' 

Using  the  definition  and  integration  by  parts,  one  obtains 


Hence, 


so  that 


A  similar  argument  on  Iq  yields 


Combining  the  above  two  equations,  we  obtain 
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.  fC2  (  3w2  3c  ,  4  \  dw 

2  =  J  (  — -T‘6cw  JT 

C1 


which  is  the  desired  formula  for  I^.  This  proves  the  lemma. 

Lemma  3.5  is  known  as  the  Picard-Fuchs  equation  for  the  integrals 
*0**1’ *2  anc*  corresPon^s  t0  Eq.  (26)  in  Cushman  and  Sanders  [10]. 


Lemma  3.6 


(3.9) 


,  2  1 

'3C  12 


3  2  ,  2 

tc  3c 

4 


Proof.  By  (3.8),  I"(c)  =  A  1 (c) (E  -  A’ (c) ) I ' (c)  where  E  denotes 
the  identity  matrix.  We  have 


A_1(c) 


2  1 

3 

1 

'  128 

1~ C 

24 

1 

£2 

1 

64  C 

8  C 

4C 

3  2 

3 

32 

32  C 

'  32  C 

2  C 

E  -  A'(c)  = 


Equation  (3.9)  follows  from  the  above  matrix  equation. 


I 


-  16  - 


Lemma  3.7  Let 


4(c) 


Then, 

(3.10) 


1  r  3  2  .  2n  1  .2. 

Q  =  8A  [TC  +  6C  Q  '  12  Q  ] 


Proof.  (3.10)  follows  from  (3.9). 

Equation  (3.10)  corresponds  to  the  Ricatti  equation  of  Cushman  and 
Sanders  [10], 


Lemma  3. 8  If  0  <  c  <  yy  »  Q  £  R  ,  then 

sH'2  *6c2Q-iq2]>0. 


Proof.  The  maximum  of  the  quadratic  form 


3  2  ,  2n  1  n2 

-  4  c  +  6c  Q  -  -lT  Q 


occurs  at  Q  =  36c  and  its  maximal  value  is 


9  r i  i , ,  2.  2 
Yj  (!  -  l44c  )c  . 


Since  0  <  c  <  yy  and  A  <  0  (Lemma  3,5),  we  have  the  desired  result. 
Lemma  3.9  If  P'(Cq)  =  0  for  some  0  <  cq  <  "fj  '  t*ien  p"(cq)  *  °* 


Proof.  By  the  definitions  of  P  and  Q,  we  have 


I'P  +  I„P'  =  I' 
0  0  2 


I"P  *  21  P’  +  1  P"  =  I” 
0  0  0  2 


I 
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:6Q  -  v2 


xoQ +  W  -  J2 

Since  P'(cQ)  =  0,  we  have  P"(cQ)  =  Iq Cc0)Q '  (cQ) / IQ (C{)) .  By  (3.6) 
and  Lemma  3.2,  Ig(c)  >  0,  I^(Cq)  <  0.  By  Lemmas  3.7  and  3.8, 

Q'  (cQ)  >  0.  This  implies  P"(Cq)  <  0. 


Lemma  3 . 10  If  P' (c^)  =  0  for  some  0  <  <  —  then 

VS  <  P(co)  <  7  • 


Proof.  By  the  first  and  last  equations  in  (3.8), 


3!0  -  16I2  -  3d*  -  12! 2 


If  P '  (Cg)  =  0,  then,  by  the  above  equation, 
3 


3  ty0 

i6-p^  =TecTM  a-4P(c0)) 


O' 


because  I£(cQ)  =  p(cQ)  1^  (cQ) .  Since  I^(c0)  <  P('C0')  “  16"  haS 

the  same  sign  as  j  -  P(Cq),  i.e.,  we  must  have  ^  <  P(cQ)  <  j  . 


Proof  of  Theorem  3.1 


It  follows  easily  from  Lemmas  3.3,  3.4,  3.9  and  3.10. 


§4.  UNIQUENESS  THEOREM  FOR  B  >_ ~  . 


In 


this  section,  we  consider  the  general  case  B  >_ j  ,  d^  =  1, 


d2  =  d^  =  d^  =  0.  From  (2.1)  the  normal  form  is 


x  =  \x  +  Bxy  +  ex 

•  1  2  2 
y  =  4  -  X  -  y 


(4.1) 


1 


-  18  - 


We  note  that  the  notations  used  in  this  section  will  be  similar  to 
those  in  §3.  This  will  not  cause  any  confusion.  In  fact,  if  we  let 
B  =  2,  then  both  notations  will  coincide.  However,  the  proofs  of 
monotonicity  in  §3  and  this  section  are  different.  The  main  reason 
to  include  the  proof  in  §3  in  this  paper  is  to  show  that,  for  B  =  2, 
one  may  be  able  to  use  complex  variable  topological  arguments  as  in 
[10,12,13]  to  prove  Theorem  3.1,  while  it  does  not  seem  likely  in  the 
present  case. 


(4.2) 


1  +  q  =  r 


Changing  time  scales,  we  obtain  the  following  normal  form 


(4.3) 


•  q  3 

x  =  xM(Xx  +  Bxy  +  ex  ) 


q,l  2  2, 

y  =  xn(-  -  x  -  y  ) 


h'!(- 


q+i  2  _  x_ 


-  X‘  y 


Using  (4.3),  we  have  along  solutions  of  (4.3) 


(4.4) 


A  =  y(x^+1X  +  ex^+3) 


Let  ey  =  X.  As  in  §3,  we  obtain,  for  periodic  solutions  of  (4.3), 
a  bifurcation  function  G(y,e,c,B)  which  for  e  =  0  is  given  by 

G(y,0,c,B)  =  (q+l)yJQ(c, B)  +  (q+3)J2 (c,B) ,  0  <  c  <  cm> 

where 
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(4.5) 


c  =  c 
m  m 


-(i)q,1(i- 7 itst)  • 

Jq(c,B)  =  f  2  xqydx  , 


J2(C 


,B)  -  f 


2  xq*2ydx  , 


•u 


_xi  c  '  18 
1+B  q+1 

X 


)  • 


and  0  <  <  a2  are  the  zeros  of  y.  We  note  that 

corresponds  to  the  heteroclinic  orbit  in  equation  (4.1) 
while  cm(B)  corresponds  to  the  fixed  point  (y ,  0)  in 
with  e  =  X  =  0.  Let 


J2(c,B) 


P^B)  =  O^B)  * 


,  0  <  c  <  c  (B) 


As  in  §3,  let 


B  B 

C1  =  al*  C2  =  a2  and 


(4.6) 


R(w) 


/  w^  w2B+2  Y5 
(  4  B+l  c  )  ’ 


0  <  c  <  c 


Note  that  R(c^)  =  R(C2)  =  0  and  0  <  c^  <  C2*  Define 


(4.7) 


VC’B^  - 


2  wnBR(w)dw,  n  =  0,1,2,, 


We  have 
(4.8) 


Mc,B) 
P(C,B)  =  l‘(c',B) 


c  =  0 

(with  e  =  X  =  0) 
equation  (4.1) 


i 


*0-  m — 
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Theorem  4.1  If  B  >  -r  and  0  <  c  <  c  ,  then 
- —  2  m 


~  (c.B)  >  0  . 

dC 


By  using  this  theorem,  it  will  be  possible  to  obtain  a  complete 
bifurcation  diagram  for  equation  (4.1)  by  following  the  methods  in 
Chow  and  Hale  [9].  We  leave  the  details  to  the  reader. 

The  following  are  needed  for  the  proof  of  Theorem  4.1. 


Lemma  4.2  For  n  =  0,1,2, 


(4.9) 


-ST  (C‘B) 


1  fc 2  wnB 


Proof.  Since 


we  have 


_  2  2B+2 

n2  w  w 

R  =  4  B+l  C  * 


2RR'  =  -1  , 


where  ”  »  ”  =  ''3/3c"  . 


Lemma  4.3  Let 


5  =  5(B)  =  lim  P (c , B) ,  n  =  n(B)  =  lim  P(c,B) 
c->0  c->-c 


Bt-I  l  = 

4  "  2  (3B+2)  <  4  T1  ‘ 


Proof.  By  (4.2),  (4.6),  (4.7)  and  (4.8),  we  have  (x  =  w  ) , 
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I2(0,B) 

Iq (0 , B) 


ii;  ^  (■  -  m 


dx 


where  = 


5  = 


fa  2 

,  .  2 

<q  (  1  - 

Jal 

V  B+l 

:  /B+l/2.  Hence 

B+l 

rv/2 

4  J 

(sin  9)' 

0 

r/2 

J 

|  (sin  0) 

n 

q+2  2, 


<)  2  c 


B+l 


2 (3B+2)  * 


Next, 


n 


31, 


, .  3c 
lim  - 

c+cm3I 


0 


3c 


(c,B) 


(c,B) 


1_ 

4 


This  completes  the  proof. 

We  are  not  able  to  find  a  differential  equation  like  (3.8)  for 
the  integrals  Iq(c,B),  I^Cc.B)  and  I2(c,B).  It  seems  this  is  possible 
when  B  is  an  integer.  Such  equations  for  B  =  1  could  be  found  in 
Carr  [8]  (Lemma  4  in  Chapter  4). 

Define 


(4.10) 


S(c 


,B)  =  [  2  w‘ 

*  r 


R(w)dw  . 
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Note  that  1^ 
check  on  some 

Lemma  4. 

(4.11) 

where  "  » "  = 


Proof. 


(c,2)  =  S(c,2).  Hence,  the  case  B  =  2  in  §3  gives 
of  the  following  formulae. 


4  The  following  equations  are  satisfied  by  1^,  S, 


‘V 

=  « 

'ro' 

B  s  • 

l' 

I' 

2 (B+2) 

0 

L  2  J 

L  2J 

"3/3c" ,  0  <  c  <  c 


m 


<f  = 


-1 


2 (3B+2) 
B+l 


and 


B+2 


c 


-1 


0 

4 


Integrating  by  parts. 
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Hence, 


The  first  equation  in  (4.11)  follows  from  the  above  equation.  The 
second  equation  in  (4.11)  is  obtained  similarly. 

Lemma  4 . 5  For  0  <  c  <  c 
-  m 


(4.  12) 


where  6  =  B/8c(B+l). 

Proof.  This  is  obtained  by  differentiating  (4.11). 


Lemma  4.6  For  0  <  c  <  c 


(4.  13) 

where 


m 


l2 

Q'  =  Q[  (1  -  4Q)  rff  +  166Q] 
2 


I’  B 

Q  =  —  and  6  =  gc(B+1) 


Proof.  This  follows  from  (4.12) 


Lemma  4.7  If  P'(Cq,B)  =  0  for  some  0  <  cQ  <  cm>  then 
P(Cq, B)  =  Q(Cq,B)  and  P"(Cq,B)  has  the  opposite  sign  as  Q'(Cq,B). 

Proof.  See  the  proof  of  Lemma  3.9. 
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Proof.  At  c  =  Cq  >  I^(cq,B)  =  P(c0,B)I^(c0,B)  since  l0  =  PIQ, 
we  obtain  (4.14)  from  the  second  equation  in  (4.11)  by  substituting 
the  above  values  of  P(Cq,B). 

Lemma  4.9  For  0  <  c  <  c 
-  m 

(4.15)  5  <  P(c,B)  <  i 

Proof.  If  not,  then  we  have  either  P(cfl,B)  _>  i  or  <_  £  for 

some  0  <  c.  <  c  .  By  Lemma  4.3,  we  may  assume  that  P'Cc^.B)  =  0. 

U  m  u 

Since  I9  >  0  and  1^  <  0,  this  contradicts  (4.14). 

Lemma  4.10  If  P'(cQ,B)  =0,  0  <  c^  <  c^,  then  P"(cQ,B)  has 

the  same  sign  as  the  following  expression  evaluated  at  (Cq,B) 

(4.16)  4V'i  *  <K?  -  KT  v2 

Proof.  By  Lemma  4.7,  P"(c0,B)  has  the  sign  opposite  to 

Q’(c0,B)  and  P(cQ,B)  =  Q(cq,B).  By  Lemma  4.9  and  (4.13),  P"(c0>B) 
has  the  sign  opposite  to  the  following  expression  evaluated  at 
(cq,B): 

u  Jfl 

(1  -  4Q)  -rf  +  !66Q 
2 
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Since  1^  <  0  and  Q(Cq,B)  <  —  ,  the  result  follows,  from  the 
definition  of  6. 


Lemma  4.11  If  P'(c0,B)  =0,  0  <  cQ  <  cm,  then  P"(cQ,B)  has 

the  same  sign  as  the  following  expression  evaluated  at  (cQ,B) 


(4.17) 


4c0r2  + 


41 J 


(^) 


Proof.  Let  P’(Cq,B)  =  0.  Using  (4.14),  we  have,  for  c  =  cQ, 

Z  4  ‘  4 

Substituting  the  above  equation  into  (4.16),  we  obtain  the  desired  re¬ 
sult. 

The  proof  of  Theorem  4.1  will  be  given  in  the  same  spirit  as  before 
(Theorem  3.1),  i.e.,  we  will  show  that  if  P'(C0>B)  =  0,  then  P"(cq,B)  >  0. 
Thus,  we  have  to  estimate  very  carefully  the  sign  of  the  expression  (4.17). 
To  begin,  define 

2  2B+2 

,  .  w  w 
r(w)  =  - - -  c  , 

2 

i.e. ,  r(w)  =  R  (w) .  Let 

/ .  N  1  -2  -2B+2 

w  "  (j)  B  and  r  =  r(w)  =  —  -  -  c  . 

Note  that  r(w)  -7  is  independent  of  c  and  3r(w)/3w  =  0.  Define 


(4.18) 


J2(c,B)  =  f  2  w2B (r (w)  -  r)  [rfw)]^  dw 
C1 
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where  c^  and  are  as  in  (4.7) 

Lemma  4.12  Let 


f(w)  =  (r  (w)  -  r) 


(  32r 

4B 

/  3r  \ 

V2  2  " 

v  3w 

w 

3w  ) 

1  - 

\  3w  / 

Then 
(4.  19) 


dw 

(f)  ^ 


Proof.  Note  that  w  is  independent  of  c  and  3r/3c  =  -1.  By 


(4.18), 

(4.20) 


J2 


2B  — 

[r (w) 


dw 


Also, 

(4.21) 


w2**(r(w)  -  r) 


(S) 


d(r3/2(w)) 


Note  that  3r/3w  has  a  simple  zero  at  w  =  w  while  the  factor  r(w)  -  r 
has  a  double  zero  at  w  =  w.  Thus,  the  following  calculations  are  valid. 
From  (4.21), 


J 


2 


(15) 


where 
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c4-22)  8M=(v  rl'4)  (r-r)+(U)  • 

'  3w  '  x  ' 

Note  that  g(w)  does  not  depend  on  c.  By  differentiating  the  above 
equation  with  respect  to  c, 


J2 


fc.  2B 
|  2  w _ 

Jr  / 


/rg(w) 


(£) 


dw 


=  -  -  I  -  r  I ' 

2  2  2 


because  of  (4.20).  Thus, 


-  —  I’  +  I '  -  r  I" 
2  2  2  2 


=  _  I  fC2  w2Bg (w) 
2  Jr  JZ 


^\W) 


2  dw 


Finally,  using  (4.9),  i.e. 


2i: 


f 


»2B  $1 


C1  [r(w)]4 


m 


dw  , 


We  obtain  (4.19)  and  that  f(w)  =  (3r/3w)  -  2g(w). 


Lemma  4.13  We  have 


(4.23) 


2rI2 


l2  " 


>  1 


Proof.  Recall  that  r  =  r(w)  and  that  3r(w)/3w  =  0.  Thus  r 
is  the  maximum  of  r(w).  From  (4.20), 
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1_ 

2 


rI2  —  ° 


Since  ^  >  0,  the  result  follows. 

Lemma  4.14  If  B  >  |  and  P’ (Cg,B)  =  0,  then  P"(cQ,B)  >  0. 

Proof.  Using  (4.23)  in  Lemma  4.11,  it  is  sufficient  to  show  that 
the  following  expression 

4c0rI2  +  ^4r  +  2c(p 1 2 

is  negative  when  it  is  evaluated  at  (c^.B).  By  Lemma  4.12,  it  is 
sufficient  to  prove  that 

h(w)  =  cQf(w)  +  (2r  +  cQ)  (9r/9w) 2  >_  0 
for  all  c^  w  c^.  Since  2B  >_  1  and  ^  0, 

2 

9  r  4B  9r 

.  2  w  9w 

9w 

2b 

=  1  -  2B  -  4w  0,  Cj  <  w  <  ^ 

By  (4.22), 

_  op  _  2 

h(w)  =  cQ(r(w)  -  r)  (1  -  2B  -  4w  )  +  2r(3r/9w) 

Since  r(w)  -  r  <_  0  and  r  ^  0,  we  have  that  h(w)  ^  0  for 
<_  w  <_  The  result  follows. 


Proof  of  Theorem  4 . 1  It  follows  from  Lemma  4.9  and  Lemma  4.14. 


§5.  BIFURCATION  FUNCTION 


Consider  a  singularity  of  type  in  the  Introduction.  This 

4 

leads  to  the  following  equation  in  R 


(5.1) 


v  =  Av  +  G(v) 


where  v  €  ]R 


A  = 


a  1 

-1  a 

0  0 

0  0 


0  0 

0  0 

B  (o 

0 


-  w 


a,0  G  1R  are  small  parameters,  w  £  R,  G  is  smooth  and  is  of  higher 
order  at  v  =  0.  If 


=  Tj  cos  0^  v2  =  -r1  sin  6^ 


v3  =  r2  cos  0  ,  v 4  =  -r2  sin  02  , 


(5.2) 


equivalent 

to  the 

following 

equations 

• 

rl 

=  «rj 

♦  R1 (0! » 02  J 

•rl*r2) 

f2 

=  0r2 

+  R2(01,02. 

®1 

=  1 

*dLW1'*2‘ 

,r1,r2) 

*2 

=  0) 

•rl’r2> 

where  each  function  is  2fr-periodic  in  0^  and  02-  Assuming  certain 
nonresonance  conditions  on  w  and  symmetric  conditions  on  G(v),  one 
can  apply  the  method  of  averaging  to  (5.2)  to  obtain  the  following 
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equations 

ri  =  rl^  ‘  aTi  '  br2^  +  •  •  * 

•  ~  2  -2 

r2  =  r2^  +  cri  +  dr2^  +  "‘ 

0j  =  1  +  . . . 

*2  =  “  +  ••• 

where  higher  order  terms  are  just  indicated  by  the  dots  and  a,b,c,d 
are  constants. 

As  in  §2,  this  leads  to  the  consideration  of  the  following  equation 
in  the  plane 

x  =  x  -  bxy2  -  x^  -  Ax  +  ex5 

(S.3) 

•  ~  2  .  3 

y  =  -y  +  cx  y  +  y 

for  x  >_  0,  y  ^  0,  where  X,e  are  small  parameters  b,  c  are  constants. 
Details  may  be  found  in  Chow  and  Hale  [9]. 

For  the  analysis  of  (5.3)  the  main  question  concerns  the  existence 
and  number  of  limit  cycles  of  (5.3)  for  (A,e)  t  (0,0)  near  the  origin. 
It  will  be  shown  that  this  is  again  related  to  the  monotonicity  of  a 
scalar  valued  function.  To  do  this,  we  must  derive  the  bifurcation 
equation. 

By  scaling  time  t  2px2^  V  and  writing 


2(c+l) 

6c-l 
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so  that  (5.3)  yields 

•  2p-l  ,  x  -2  3. 

x  =  2px  r  y(-x  -  bxy  -  x  ) 

(5.4)  +  x^’Vc^x  +  e^5)  , 

y  =  2px2p_1y(-y  +  cx2y  +  y3)  , 

where 

-  ^+1  T  -  *  +  b 

p  ~  Be- 1  ’  =  6c-l  ' 

For  e1  =  =  0,  the  equation  (5.4)  has  a  first  integral.  This 

can  be  used  to  obtain  the  bifurcation  function  for  periodic  orbits. 

We  have  not  been  able  to  analyze  completely  the  bifurcation  function 
for  this  case.  There  seem  to  be  regions  in  the  (b,c) -parameter  space 
where  there  is  only  one  periodic  orbit  and  other  regions  where  there 
are  more  than  one  periodic  orbit.  Therefore,  to  illustrate  non¬ 
uniqueness,  we  assume  in  the  remainder  of  this  section  that  x  =  1, 
p  >  0.  These  relations  are  equivalent  to  b  =  2/ (c  -  1) ,  c  >  1 .  Notice 
that  this  implies  that  be  >1. 


(5.5) 


3H 

8X 


=  0,  we  have  the  following  Hamiltonian  system: 
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Fig.  5.1 


dH 

dt 


2P-1  .  5,  3H 

x  y(XjX  +  ejX  )  — 


^  5,  2p-l  d 

-  (XjX  +  EjX  )  x  v  ^ 


Integrating  by  parts  from  time  t^  to  time  t^, 

H(x(t2),y(t2))  -  HCxftp  ,y(t1)) 


(X^x  +  e1x5)x2p  1 


x(t.) 

x(tj)  2 


2pX1x2p~1  +  Ej  (2p  +  4)x 


2p+3 


dx 


=  L(x,y) 


x(t2) 


xitp 


M(x,y)dx,  say. 


Let  t1  =  0.  Assume  (x(0),y(0))  is  on  the  line  l  defined  in 
(5.10).  Let  t2  >  0  (t”  <  0)  be  the  first  forward  (backward)  time  that 

the  periodic  orbit  (x(t),y(t))  meets  the  line  l.  Hence, 


i 
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H(x(t*),y(t*))  -  H(x(0) ,y(o)) 
tt  fx(t|) 

=  L(x,y+)  +  M(x,y+)dx 

-  0  Jx(0) 

+  oc|x1  +  £l|2), 

where  y+  is  given  by  (5.7).  Thus,  periodicity  of  (x(t),y(t)),  i.e., 
HfxCt^.yCtp)  =  H(x(t~),y(t')),  is  equivalent  to 

fx(t*)  fx (t ”) 

(5.10)  M(x,y  )dx  =  M(x,y_)dx 

■’x(O)  x(0) 

*  0(|»,  *  ej2)  . 

By  (5.7)  and, by  using  the  correct  signs,  (5.10)  is  equivalent  to  the 
following: 

0  =  |  2  ATx)  |^2pXjX2^_1  +  Ej(2p  +  4)x2p+3j  dx 
Jxl 

+  0(|Xj  +  ej2) 

where  x^  =  x(0)  and  =  x(t*)  =  xCt^).  Let  =  eiyi’  By 
dividing  the  above  equation  by  we  obtain  the  bifurcation  equation 

0  =  G(y1,c1,c) 

=  p  ATxT  [2pu1x2p_1  ♦  (2p  +  4)x2p+3]  dx 
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By  (5.8)  and  (5.9), 


(5.11) 


where 

(5.12) 


G(p1,0,c) 


(p  +  2)w 


£  + 

2 


‘] 


R(w)dw 


R(w)  =  [w^(pw  -  l)2  -  4c] 


Define 


Jn(c) 


^ -l+2n 

w  R(wjdw, 


n 


0,1  . 


Consider  the  function 

J,(c) 

(5.13)  Y(c)  = 


We  summarize  the  above  discussion  into  the  following. 


Theorem  5 . 1  If  (A^,e^)  is  sufficiently  small,  then  every  periodic 

orbit  of  (5.5)  must  intersect  the  line  *.  defined  by  (5.10)  with 

0  <  x  <  1/2+p.  Furthermore,  there  exists  a  smooth  function  yjfe^c) 

such  that  a  necessary  and  sufficient  condition  for  (x(t),y(t))  to  be 

a  periodic  solution  of  (5.5)  with  (x(0),y(0))  £  l,  0  <  x2(0)  <  1/2+p, 

H(x(0),y(0))  =  c,  is  that  =  ^(e^c),  where  A^  =  .  Further- 

more,  2 

wj(0,c)  =  -Y(c)E-£  . 

To  discuss  the  behavior  of  the  function  Y (c)  as  a  function 
of  c,  observe  that  the  range  of  values  of  c  is  [0,c*]  where 


c*  >  0  and 
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(5.14)  W]i(0)  =  0,  w2«J)  =  i 

(5.15)  wi(c*)  =  w2(c*)  = 

Using  (5.14),  one  can  easily  integrate  Jj(0),Jq(0)  to  obtain 

Y(0)  =  - Eli - 

1  J  P (P  +  2)  (p  +  3) 

Using  (5.15),  one  observes  that 

Y  (c*) - 

(P*2)Z 

and 

(5.16)  Y (0)  -  Y (c*)  =  - ^ -  >  0 

p(p+2)  (p+3) 

Numerical  evidence  indicates  that  Y(c)  is  monotone  and  thus  periodic 
orbits  are  unique.  (See  Fig.  5.2.) 

It  would  be  interesting  to  obtain  the  bifurcation  function  for 
t  1  and  also  to  determine  numerically  or  analytically  if  there  are 
regions  in  parameter  space  where  there  is  nonuniqueness  of  periodic  orbits. 
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APPENDIX  A 

If  B  =  2,  then  the  normal  form  is 


(Al) 


•  3  2 

x  =  Xx  +  2xy  +  efdjX  +  c^xy  J 

*  1  2  2  2  ,  3, 

y  =  j  -  y  -x  +  e(d3x  y  +  d4y  ) 


where  d^,  d^,  d^  and  d^  are  constants.  In  §2  and  §3,  we  assumed  that 
dj  =  1,  d2  =  d3  =  d^  =  0.  As  in  [8]  or  [9],  we  will  obtain  a  bifurca¬ 
tion  function  G(y,e,c)  which  for  e  =  0  is  given  by 


(A2) 


G(y,0,c)  =  y 


r&2 

ydx  +  K 
ai  2 


t  a 


^  x^ydx  +  K 
a-l  0 


/•a. 


ydx 


where  X 


(A3) 


=  yJQ  ♦  K2J2  ♦  K0Jq 


ey,  JQ,  J2  are  defined  by  (3.2)  and  (3.3)  and 


IC2  .  JWj  -  d4)  -  Cd2  -  d3) 


2  3 

K0  =  T  +  4  d4 


To  show  (A2) ,  we  note  that  an  orbit  r  of  (Al)  is  periodic  with 
minimal  period  T  >  0  if  and  only  if 


where 


0  = 


H  dt, 

r 


and  (x,y)  is  the  periodic  solution  T  of  (Al).  By  using  X  =  ye 
one  obtains  the  following 
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2G(y,0,c) 


•  3  2 

y(yx  +  d^  +  d2xy  )dt 


-fi 

}0 


x(d3x2y  +  d4y3)dt 


where  (x(t),y(t))  is  the  solution  of 


x  =  2xy 


•  1  2  2 
X  =  4  -  x  -  y 


with  energy  H(x(t),y(t))  =  c.  Since  all  the  terms  in  (A4)  are  given  by 
JQ  and  (see  (3.2)  and  (3.3))  except  the  following  expressions 


1  fT  3..,  fc 
7  y  xdt  = 
z  J0  £ 


2  y3dx 


1  fT  +  1  fa2 

2  xy  ydt  =  -  -  y  dx  , 
J0  J 


we  will  show  that 


can  be  expressed  in  terms  of  JQ  and  J2- 


By  (A5), 


2xy  ^  -  x2  -  y2 

y  dx  4  X  V  • 


1 
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Integrating  by  parts. 


2 

y  dx 


■  f2  -  V  £  a. 


3  Ca2  fl 

5  .  yt4  ■ 


2  2 
x  -  y  )dx 


-  7  Jo  -  3J2  • 


This  gives  (A2)  and  (A3) . 

If  K2  /  0,  then  by  Theorem  3.1  we  have  a  unique  asymptotically 
stable  limit  cycle  of  (Al)  for  appropriate  values  of  (A,e)  ^  (0,0). 
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APPENDIX  B 

If  B  i  2,  the  normal  form  is  given  by 


(Cl) 


•  3  2 

x  =  Ax  +  Bxy  +  e(d^x  4  d^xy  ) 

•  1  2  2  2  ,3. 

y  =  _  .  x  -  y  +  e(d3x  y  +  d4y  ) 


where  d^,  d^,  d^  and  d^  are  constants.  Changing  time  scales,  we 
write  (Cl)  as 


(C2) 


Let 


x  =  x^[Ax  +  Bxy  +  efd^x^  +  d^xy2)] 

•  qrl  2  2  ,  ,  2  3., 

y  =  xM[^  -  X  -  y  +  c(d3x  y  +  d4y  )]  . 


H  =  I  .  xq+1y2 .  I 

24  y  1+B  j  * 


where  1+q  =  2/B.  Using  (C2),  it  follows  that 

dH  •  q+l  •  3  ,  2.,  q 

^  =  yxM  A  +  ey(dxx  +  d2xy  )xM 


•  , ,  2  ,  3,  q 

-  ex(d3x  y  +  d4y  )xM 


As  in  Appendix  A,  the  bifurcation  function  G(p,e,c,B),  ep  =  X,  which 
for  e  =  0  is  given  by 

G(p,0,c,B) 

"  ^+1^J0  +  K2J2  +  Vo  * 


where 


<2  =  (q+3)d1  -  (q+l)d2  4  d3  -  3d4  , 


K0  =  (d2(q+D  4  3d4)/4  , 

and  and  J2  are  defined  in  §4.  This  shows  that  there  is  no  loss 

of  generality  by  working  with  equation  (4.1). 


i 
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